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The thermoelectric power, including the Nernst and Seebeck effects, in graphene nanoribbon is studied. By 
using the non-equilibrium Green function combining with the tight-binding Hamiltonian, the Nernst and See- 
beck coefficients are obtained. Due to the electron-hole symmetry, the Nernst coefficient is an even function of 
the Fermi energy while the Seebeck coefficient is an odd function regardless of the magnetic field. In the pres- 
ence of a strong magnetic field, the Nernst and Seebeck coefficients are almost independent of the chirality and 
width of the nanoribbon, and they show peaks when the Fermi energy crosses the Landau levels. The height of 
rc-th (excluding n = 0) peak is [ln2/|«|] for the Nernst effect and is ln2/« for the Seebeck effect. For the zeroth 
peak, it is abnormal with height [2 In 2] for the Nernst effect and the peak disappears for the Seebeck effect. 
When the magnetic field is turned off, however, the Nernst effect is absent and only Seebeck effect exists. In 
this case, the Seebeck coefficient strongly depends on the chirality of the nanoribbon. The peaks are equidistant 
for the nanoribbons with zigzag edge but are irregularly distributed for the armchair edge. In particular, for the 
insulating armchair ribbon, the Seebeck coefficient can be very large near the Dirac point. When the magnetic 
field varies from zero to large values, the differences among the Seebeck coefficients for different chiral ribbons 
gradually vanish and the nonzero value of Nernst coefficient appears first near the Dirac point then gradually 
extents to the whole energy region. 

PACS numbers: 72.15.Jf, 73.23.-b, 73.43.-f, 81.05.Uw 



I. INTRODUCTION 

As a single atomic layer extracted from graphite, graphene 
has been successfully fabricated experimentally i 2 - Due to its 
peculiar topological structure, the graphene exhibits peculiar 
properties. 3 - For the graphene sheet, the conduction and va- 
lence band in graphene intersect at Dirac points, the corners 
of the hexagonal first Brillouin zone. Around the Dirac points 
graphene has a unique band structure and its quasi particles 
satisfy the massless Dirac equation where the speed of light is 
replaced by the Fermi velocity of graphene (\'f m I0 6 m/sec). 
Experimentally, by varying the gate voltage, the charge car- 
riers of graphene can be easily tuned, globally^ or locally^. 
As a result the Fermi level can be above or below the Dirac 
points, which is viewed as electron-like or hole-like system. 
Along the different crystal direction in honeycomb lattice, 
the band structure 5 and the transport properties are different. 
For the graphene ribbon with the zigzag edge, a special edge 
state exists^ While for the graphene ribbon with armchair 
edge, it is metallic when the transverse layer N = 3M - 1 
with integer M and insulator otherwise^ When the perpen- 
dicular magnetic field is strong enough to form Landau lev- 
els (LLs), these differences due to different chirality at the 
zero magnetic field disappear. In addition, both theoretically 7 
and experimentally^ the Hall conductance was found to be 
the half-integer in the values g(n + 1 /2)e 2 /h with degeneracy 
g = 4, indicating that the quantization condition is shifted by a 
half-integer compared with the usual integer quantum Hall ef- 
fect. It is a direct manifestation of the unique electronic struc- 
ture of graphene. 

The thermoelectric power (TEP), or the thermal gradient in- 
duced current (or bias with an open boundary), results from a 
balance of electric and thermal forces acting on the charge car- 
riers. In general, we consider two thermoelectric powers, the 



Nernst effect which is the transverse TEP induced by a lon- 
gitudinal thermal gradient in a perpendicular magnetic field 
and the Seebeck effect which is the thermal gradient induced 
bias in a two probe system. TEP is of great importance in un- 
derstanding electronic transport because it is more sensitive 
to the details of the density of states 8 , and the particle-hole 
asymmetry 9 than the conductance. In the early days, because 
of the experimental difficulty (particularly in low-dimensional 
systems or nano-devices), the Nernst effect and Seebeck effect 
are often neglected. Instead, one usually measures the Hall 
effect and the resistivity. Now, with the development of the 
micro-fabrication technology and the low-temperature mea- 
surement technology, the thermoelectric measurement in low- 
dimensional samples has been feasible^ Recently, the Nernst 
effect and Seebeck effect have been widely observed and 
experimentally investigated in many systems, including the 
high-Tc superconductivity, 1 - ferromagnets^ semimetallicr^ 
Pyrochlore Molybdates, 14 Bismuch, 15 single walled carbon 
nanotube^ etc. For the graphene, the study of thermoelectric 
properties can elucidate details of the electronic structure of 
the ambipolar nature that cannot be realized by probing con- 
ductance alone. Very recently, using a microfabricated heater 
and thermometer electrodes, the conductance and the diffu- 
sive TEP of graphene are simultaneously measured by Zuev 
et.alM- and Wei et.aJ.J^ Zuev et.al. found electrons and holes 
contribute to Seebeck effect in opposite ways. At high temper- 
atures direct measurement of Seebeck coefficient Sc can be 
compared with that calculated from the Mott relation.— Fur- 
thermore, divergence of Sc and the large Nernst signal were 
found near the charge neutral point (i.e. the Dirac point). 18 
Also, at low temperatures, depending on £>, TEP is oscillat- 
ing. The temperature suppresses the oscillation and enhances 
the magnitude of TEP. 

Up to now, some theoretical investigations have been car- 
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ried out on the thermal response in the graphene. The elec- 
tronic transport coefficients including thermopower was semi- 
classically treated and only classical Hall effect (low field) in 
graphene was studied, 20 In addition, the Nernst coefficient was 
studied only in weak magnetic field. It was found to be strong 
and positive near Dirac point.— For a strong magnetic field in 
the quantum Hall regime, the Seebeck coefficient was studied 
and was focused on its dependence of the field orientation. 2 ' 
In all these works, the quantum Nernst effect is absent be- 
cause of the calculational subtleties in the presence of the 
strong magnetic field. For the normal two dimensional elec- 
tron gas characterized by a parabolic dispersion, the Nernst 
effect has been studied^— Of these works, two alternative 
boundary conditions were considered in calculating the ther- 
mal response functions. One is the adiabatic boundary con- 
dition that the temperatures in the upper and lower edge are 
fixed. In this case the Nernst coefficient is similar to the See- 
beck coefficient. 22 The other one is the non-adiabatic bound- 
ary condition on the upper and lower edges, in which the 
edge currents are in contact with two heat baths with differ- 
ent temperatures, 24 The Nernst coefficient is different from the 
Seebeck coefficient. It is the purpose of our work to focus on 
the quantum Nernst effect in the graphene nanoribbon with 
the adiabatic boundary condition. 



In this paper, we carry out a theoretical study of the Nernst 
effect in a crossed graphene nanoribbon and the Seebeck ef- 
fect in a single graphene nanoribbon in the strong perpendic- 
ular magnetic field, zero magnetic field, and weak magnetic 
field. By using the tight binding model and the nonequi- 
librium Green function method, the transmission coefficient 
and consequently the Nernst and Seebeck coefficients are ob- 
tained. In a strong perpendicular magnetic field B, high de- 
generated LLs are formed, and the edge states dominate the 
transport processes, so the Nernst (Seebeck) coefficients are 
almost the same along the different chiral directions. We find 
that the Nernst coefficient Nc and the Seebeck coefficient Sc 
show peaks when the Fermi energy Ef passes the LLs. At 
Ef = 0, because the zeroth LL is shared by electron-like and 
hole-like Landau states, Nc which is an even function of Ep 
has the highest peak while Sc which is an odd function of Ef 
vanishes. On the other hand, at zero B, there is no Lorentz 
force to bend the trajectories of the thermally diffusing carri- 
ers, so Nernst effect is absent. In this case, the Seebeck coef- 
ficient S c is strongly dependent on the chirality of graphene 
ribbon. In particular, for the insulating armchair ribbon, Sc 
can be very large near the Dirac point. At last, the crossover 
behavior of the thermoelectric power from the zero magnetic 
field to the strong magnetic field is also studied. 



The rest of the paper is organized as follows. In Section II, 
the models for crossed graphene ribbon or single graphene rib- 
bon are introduced. The formalisms for calculating the Nernst 
and Seebeck coefficient are then derived. Section III gives nu- 
merical results along with discussions. Finally, a brief sum- 
mary is presented in Section IV. 



II. MODEL AND FORMALISM 

We consider two graphene systems: a four terminal crossed 
graphene nanoribbon and a two terminal graphene nanoribbon 
as shown in the left and right insets of Fig. 1(b). Here we con- 
sider ballistic two dimensional electron gas in which the mean 
free path and the phase coherent length are greater than the de- 
vice size. In the experiment, we can use the smaller sample to 
reduce the device size, and use the lower temperature or the 
higher magnetic field to enhance the phase coherent length. In 
the tight-binding representation, the Hamiltonian operator can 
be written in the following form ; 3 ' 25 ' 26 



= ^ eifl.'fli - 2^ fe'^flj'flj, 
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where i = (i x , i y ) is the index of the discrete honeycomb lattice 
site which is arranged as in inset of Fig. lb, and a\ and a\ are 
the annihilation and creation operators at the site i. e, is the 
on-site energy (i.e. the energy of the Dirac point) which can 
be controlled experimentally by the gate voltage, here we set 
e, = as an energy zero point. The second term in Eq.(2) is the 
hopping term with the hopping energy t . When the graphene 
ribbon is under a uniform perpendicular magnetic field B z = 
B, a phase 4>ij is added in the hopping term, and 0y = J. J A ■ 

dl/00 with the vector potential A = {-By, 0,0) and the flux 
quanta <po = H/e. 

With this ballistic system, the current flowing to the p-th 
graphene lead can be calculated from the Landauer-Biittiker 
formula^ 

J P = jj] f ^{T pq (E)(f p (E)-f q (E))]. (2) 

i 

where p, q = 1,2,3,4 for the four terminal system or p, q = 
1 , 2 for the two terminal system, and T pq is the transmission 
coefficient from terminal-q to terminal-p. 

In Eq.©, the transmission coefficient T pq can be calculated 
from T pq (E) = J r[T P G'T q G a ], where the line-width function 
T p (E) = /(I;;-!);'). The Green's function G'\E) = [G a (E)] f = 
\EI-Hq-Y> p £p(£)L' where Hq is Hamiltonian matrix of the 
central region and / is the unit matrix with the same dimen- 
sion as that of Hq, and V p is the retarded self-energy function 
from the lead-p. The self-energy function can be obtained 
from T. r p (E) = H c , p g r p (E)H PtC , where H cp (H PtC ) is the cou- 
pling from central region (lead-p) to lead-p (central region) 
and g' p {E) is the surface retarded Green's function of semi- 
infinite lead-p which can be calculated using transfer matrix 
method3/ p (£') in Eq.(2) is the Fermi distribution function, it 
is also a function of the Fermi energy Ef and temperature T , 
and can be written as 



f p {E,E p F ,T p ) 



1 



„(E-E p F )/k B T„ + 1 



(3) 



where E P F — Ef + eV p with e the electron charge and V p 
is the external bias. In the four terminal system, the ther- 
mal gradient AT is added between the longitudinal terminal- 
1 and terminal-3, and 71 = T + 0.5A7", T 3 = T - 0.5A7", 
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V\ — V3 = 0. Due to the Lorentz force, the longitudi- 
nal thermal gradient induces a transverse current J 2 ^ m the 
closed boundary condition or a transverse bias V 2 ,4 in the open 
boundary condition in the terminal-2 and terminal-4. Here we 
consider the open boundary (72 = V4 = 0) and calculate the 
balanced bias V 2 a- While in the two terminal system, both 
original thermal gradient AT and induced balanced bias are 
considered in the longitudinal terminal- 1 and terminal-2, and 
we have T\ - T + 0.5 AT and T 2 = T - 0.5 AT. Assuming 
small thermal gradient and consequently the small induced ex- 
ternal bias, the Fermi distribution function in Eq.© can be 
expanded linearly around the Fermi energy Ef and the tem- 
perature T, 



f P (E,E P F ,T p ) = f + eV p 



= /o+/o(/o-D 



AT, 



df 



V P =Q,T„=T 

eV, 



p dT B 



V„=0,T„=T 



k R T 



+ {E-E F ) 



AT,, 



k B T 2 



where fo = |e (£ e f)/i<bT + j s me Fermi distribution in the 
zero bias and zero thermal gradient. Then for the four terminal 
system, the current J 2 of the terminal-2 can be rewritten as: 



h = 



^ J" dEf (f -l)T 21 (E) 
^JdE Mf - l)T 23 (E) 
y JdE /„(/,) - l)T u (E) 



2k B T 2 k B T 
-AT qV 2 



2k B T 2 k B T 



v 2 -v 4 



k B T 



(5) 



Similarly, the expression for the current J 4 of the terminal-4 
can also be obtained. Using the open boundary condition with 
J2 = J a — and considering the system symmetry (T 2 \ = T43, 
T23 = T41 and T 2 4 — T42), the Nernst coefficient Nc in the 
four terminal system is: 



N c 



V2-V4 
AT 

1 JdE(E- E F )(T 2l - r 23 )/o(/o - 1) 
eT 5dE{T2x+T23+2T24)f a {fo-l) 
In the two terminal system, the current J\ = -J 2 is 



(6) 



./ 



-/ 



dEf (f - l)T 21 (E) 



2k B T 2 

Let J 1 = 0, we have Seebeck coefficient S c 
V1-V2 



k B T 



S c = - 



AT 



1 f dE (E - £»r 2 i (£)/()(! - fo) 
JdET 2l (E)Ml-f ) 



(7) 



in. NUMERICAL RESULTS AND DISCUSSION 

In the numerical calculations, we set the carbon-carbon dis- 
tance a = 0.142«m and the hopping energy t = 2.75eV as 
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FIG. 1: (Color online) Nernst coefficient Nc (a) in the four termi- 
nal system and Seebeck coefficient 5 c (b) in two terminal system vs 
Fermi Energy E F with the strong magnetic field BSq = O.OOS(po/n 
and ribbon width N = 80. Different curves are for different tempera- 
tures k B T. The four terminal system and the two terminal system are 
shown in left and right inset in panel (b), respectively. 



in a real graphene sample. 4 ' 5 Throughout this paper the en- 
ergy is measured in the unit of t. The magnetic field B is 
expressed in terms of magnetic flux BSo in the unit of <f>o/n 
where So = § V3a 2 is the area of a honeycomb unit cell and 
0o = h/e is the flux quanta. If we set BSo - 0.00 l<f>o/n, the 
real magnetic field is around AT . The width of the graphene 
ribbon is described by an integer N, and the corresponding 
real width is 3Na for zigzag edge nanoribbon and y/3Na for 
the armchair edge nanoribbon. In the schematic setup- II and 
setup-Ill in the inset of Fig. lb, N - 2. In the presence of the 
strong perpendicular magnetic field, since transport properties 
are independent of the chirality, we choose the setup- II shown 
in the left inset of Fig. 1(b) to study the Nernst effect and the 
setup-111 shown in right inset of Fig. 1(b) to study the Seebeck 
effect. On the other hand, when the magnetic field is zero, 
the Seebeck effect strongly depends on the edge chirality, so 
we will study both zigzag and armchair edge nanoribbons, re- 
spectively. 



A. the strong perpendicular magnetic field case 

Firstly, we study the system with strong perpendicular mag- 
netic field. Fig. 1 shows the Nernst coefficient Nc and Seebeck 
coefficient S c versus Fermi energy E F for different tempera- 
tures T = 0.001 f, 0.003r, 0.006f and O.Olf. Considering the 
ambipolar nature of the graphene and the electron-hole sym- 
metry, the Nernst coefficient Nc is an even function of Ef 
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[Nc(Ef) = Nc(-Ef)], because both the energy E—Ef and the 
direction of the particle movement (or Ti\ - T23) reverse their 
signs under the electron-hole transformation. From Fig. 1(a), 
we see that the Nernst coefficient Nc show peaks when Ef 

passes the LLs E„ = sign{n) ^j2efiv 2 F \n\B and show valleys 
between adjacent LLs. With the increase of the temperature, 
the peak heights roughly remain unchanged, but the valleys 
rise. For convenience, the peaks are numbered and the peak at 
Ef — is denoted as the zeroth peak. In the low temperature 
limits, for the n-th peak with n + 0, the height is [In 2/|n|], and 
the zeroth peak height is [21n2]. In Fig. 2(a) we plot inverse 
of the peak heights versus the peak number n (see the crossed 
circle symbols) at the low temperature T = O.OOlf. It satis- 
fies the relation [|«|/ ln2]. For comparison, the inverse of the 
peak's height for the conventional metal is also plotted (see 
dotted pentagram symbol), which is [(« + 1 /2)/ ln2]. 

In Fig. 1(b), we plot the Seebeck coefficient Sc versus Ep 
at different temperatures T. Similarly, the Seebeck coefficient 
Sc display peaks when Ef passes the LLs and show valleys 
between adjacent LLs. However, Sc shows two essential dif- 
ferences from the Nernst effect: First, Sc is an odd function 
of Ef, which means that contributions to Sc from electrons 
and holes differ by a sign due to the electron-hole symmetry. 
So the Seebeck coefficient Sc is negative for Ef < 0. Second, 
when Ef is on the zero-th LL, S c is zero instead of the highest 
peak in the curve of Nc-Ef. This is because the zero-th LL 
with the fourfold degeneracy is shared equally by electrons 
and holes and the electrons and holes give the opposite con- 
tributions to Sc- The inverse of the peak height of Seebeck 
coefficient at the low temperature (ksT = O.OOlf) is plotted 
in Fig. 2(b). It is found that in graphene, the pseudospin re- 
lated Berry phase 2 introduces an additional phase shift in the 
magneto-oscillation of TEP. As a result of this phase shift, the 
inverse of peak height is oc « (see the crossed circle symbol 
in Fig. 2(b)). While in the conventional metal or semiconduc- 
tor with massive carriers, there is no pseudospin related berry 
phase, the inverse of peak height is oc « + ± (see dotted pentacle 
symbol in Fig. 2(b)). 

Next, we study the temperature effect. Since TEP (Nernst 
effect or Seebeck effect) represents the entropy transported 
per unit charge, both Nernst coefficient and Seebeck coeffi- 
cient increase with the increasing temperature which are ex- 
hibited in Fig. 1(a) and (b). To take a closer look in Fig. 3, 
we plot the zeroth and first peak for the temperature range 
6 [0.001f,0.029f] in the step of 0.002f. The temperature ef- 
fect of Sc is similar to that of Nc, so we only show the Nernst 
coefficient Sc in Fig. 3(a). At low temperatures, with the in- 
crease of the temperature kgT, the peak height and position 
do not vary much, but the peak half-width is broadened pro- 
portional to ksT, so the valley between the LLs rises. When 
the temperature kgT exceeds the spacing of nearest LLs, the 
Nernst and Seebeck coefficients Nc and S c are enhanced in 
the whole range of energies including both the peak and val- 
ley because of the overlap of the neighboring peaks. In ad- 
dition, except for the zeroth peak, the peak positions for all 
other peaks shift towards the zeroth peak. 

Now we study the disorder effect on the Nernst and See- 
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back effect. To consider the effect of disorder, random on-site 
potentials (5e, in the center region are added with a uniform 
distribution [-W/2, W/2] with disorder strength W. The data 
is obtained by averaging over up to 1200 disorder configu- 
rations. It is known that when the magnetic field is absent, 
the Seebeck effect is strongly affected by the disorder, and the 
peaks are suppressed even in the small disorder. On the other 
hand, in the presence of the strong magnetic field, the Seebeck 
effect and Nernst effect are robust to the disorder, because of 
the existence of the quantized Landau level. The bigger the 
sample is (or the stronger the magnetic field is), the more ro- 
bust the Nernst effect and Seebeck effect. Similar to Fig. 3(a), 
in Fig. 3(b) we plot the zeroth and first peak at fixed tempera- 
ture kgT = O.Olf with different disorder strengths. Here sam- 
ple size (N = 40) is smaller than that in Fig. 1(a) (in which 
N = 80). With the smaller sample size, the zeroth universal 
values of peak height 2 In 2 can still remain until disorder W 
is larger than If. For the first peak, the universal values of 
height In 2/|n| remains at W — 0.3f and washes out at stronger 
disorder. It means that the Nernst peak corresponding to the 
lower Landauer level can resist stronger disorders. In fact, this 
effect of disorder has been studied for the thermal response to 
the charge current^ or to the spin current^ So, in the follow- 
ing, we will focus only on the clean system. 



B. the case of zero magnetic field 
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In this subsection, we study the TEP at zero magnetic field. 
Because there is no Lorentz force to bend the trajectories of 
the thermally diffusing carriers, the Nernst effect is absent and 
Nc = 0. At B = 0, the Seebeck coefficient S c is strongly de- 
pendent on the chirality of graphene ribbon. In addition, for 
the armchair edge ribbon, it is metallic when N = 3M - 1 (M 
is an integer) and insulator otherwise; 6 The Seebeck coeffi- 
cient S c has essential difference for the metallic and insulator 
armchair ribbons. In the following we consider three different 
systems: (1) zigzag edge ribbon with width N = 40 (sketched 
in inset of Fig. 4(a)), (2) metallic armchair edge ribbon with 
width N - 41 (sketched in inset of Fig.4(b)), and (3) insulat- 
ing armchair edge ribbon with width N = 40 (sketched in inset 
of Fig.4(b)). Fig.4(a), (b), and (c) show the Seebeck coeffi- 
cient S c versus Ef for the above three systems, respectively. 
For the convenience of discussion, we also plot correspond- 
ing transmission coefficient T — Ti\ — T\o versus Ef in each 
panel. We can see that S c is an odd function of Ef and S c in- 
creases when the temperature increases. In addition, S c peaks 
when Fermi energy crosses the discrete transverse channels 
where quantized transmission coefficient jumps from one step 
to another. These properties are similar for the above three 
cases. 

But there are also many essential different behaviors. (1). 
For the zigzag edge ribbon, the transverse channels are 
equidistance with the energy interval A = |f|7r/(2A0 in the 
conduction band or the valence band (except that the interval 
from the first transmission channel in the conduction band to 
the first transmission channel in the valence band is 3A). So 
peaks of S c are uniformly distributed over energies and the 



FIG. 4: (Color online) Seebeck coefficient Sc vs. Fermi energy E F 
for the different temperatures k B T at zero magnetic field. Panel (a) 
is for the zigzag ribbon as sketched in the inset of panel (a), with the 
width N = 40. Panel (b) and (c) are for the armchair ribbon sketched 
in inset of panel (b) with the width N = 41 (b) and N = 40 (c). The 
gray solid curves in panels (a), (b), and (c) are the corresponding 
transmission coefficients T. 
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FIG. 5: (Color online) The Seebeck coefficient of Fig.4(c) with the 
Fermi Energy interval [-0.03, 0.03]. 
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80. 



peak height of Sc satisfies [ln2/2n] where « is the peak num- 
ber (see Fig.4(a)). (2). In metallic armchair edge ribbon, how- 
ever, the transverse channel and consequently the peaks of S c 
are irregularly distributed. The peak height of S c is closely 
related to the transmission coefficient T = T%\ - T\i and it 
can be expressed as 2ATln2/(2T + AT) at low temperatures, 
where AT is the change of T when Ef scans over the certain 
transverse channel. With increasing of the temperature, some 
of peaks that are very close to each other merge together so 
that both peak height and position are irregular (see Fig.4(b)). 
(3). Finally, for the insulating armchair edge ribbon, except 
for the irregularly distributed peaks for \Ep\ > A, the Seebeck 
coefficient Sc is very large for Ef near the Dirac point (0) at 
low temperatures. Fig. 5 magnifies the curves of Sc-Ef near 
the Dirac point. At low temperatures, S c can be very large 
when Ef approaches the Dirac point. For example Sc can 
reach about 10 at T — 0.0022?. At the Dirac point the sign of 
S c changes abruptly. This is because near the Dirac point the 
transmission coefficient T\i is zero and the carriers can't be 
transmitted. In order to balance the thermal forces acting on 
the charge carriers, we have to add a very large bias leading 
to a very large Seebeck coefficient near the Dirac point at low 
temperatures. When temperature increase such that k^T is 
greater than the gap of the insulating armchair edge ribbon S c 
decreases gradually. We emphasize that if the armchair edge 
ribbon is narrow enough (such as W ~ \0nm as in our calcu- 
lation), S c ~ 10 at the temperature T = 0.0022f/jfc B ~ 60K. 
This very large S c can be observed in the present technology. 



the crossover from zero magnetic field to high magnetic 
field 



In this subsection, we study the Nernst and Seebeck ef- 
fect when the magnetic field varies from zero to finite values 
(strong magnetic field). At zero magnetic field, the Nernst co- 
efficient Nc is zero and the Seebeck effect S c is dependent 
on the chirality of graphene ribbon. At high magnetic fields, 
however, both Nc and S c are independent of the ribbon chiral- 
ity. What happens with the magnetic field in the intermediate 
range? 

First, we study the Nernst effect, in which two different se- 
tups (the setup-61 and setup-611) sketched in the top of Fig. 6 
are considered. In Fig. 6 we plot the Nernst coefficient Nc ver- 
sus Ef at different temperatures and magnetic fields. From 
Fig. 6(a) to (c), the magnetic field increases from weak to 
strong enough to form edge state. At the weak magnetic field 
(such as BS[)/(po = 0.0005/7r), the Nernst coefficient Nc peaks 
sharply near the Dirac point at low temperatures. Because 
on two sides of the Dirac point, the carriers are electron-like 
and hole-like and they are shifted to the opposite direction un- 
der the weak magnetic field, the Nernst effect is largest at the 
Dirac point. 

In particular, in the setup-61, the Nernst coefficient Nc is 
very large at the Dirac point, which is much larger than that 
in setup-611 and in the case of high magnetic field. Because 
for the setup-61, the longitudinal leads (lead-1 and lead-3) are 
metallic with a large transmission coefficient but the trans- 
verse leads (lead-2 and lead-4) are almost insulator near the 
Dirac point. As a result, we have to add a much larger bias 
to balance the thermal current so that the Nernst coefficient 
Nc is very large in the setup-61 at the low magnetic field (see 
Fig.6(al) and (bl)). With increasing of B, LLs are formed one 
by one. The zeroth LL located at the Dirac point is formed 
first (at about B5o/0o = 0.001 5 In, no shown), then is the first 
LL, the second and so on. For example, In Fig. 6(a), no LL 
is formed while in Fig. 6(b), the zeroth, first and second LL 
are formed. As soon as LLs are formed, the Nernst coefficient 
Nc will satisfy the relation that its peak heights are equal to 
In 2/|n| (or 2 In 2 for n = 0). From Fig. 6(c), we can see that 
as BSo/<po = 0.005 /7r, electrons (or holes) with Fermi energy 
\Ef \ < 0.3f all belong to robust edge states. In this case, the 
Nernst coefficient Nc are almost the same for the setup-61 and 
setup-611. 

For the Seebeck effect, armchair edge ribbon can either be 
metal or insulator, we also consider three different systems as 
in the case of zero magnetic field. In Fig. 7 we plot the See- 
beck coefficient Sc versus Ef at different temperatures and 
magnetic fields for three different systems. The first column 
is for the zigzag edge ribbon with width N = 80, the second 
column is for metallic armchair edge ribbon with N = 80, 
and the third column is for insulating armchair edge ribbon 
with N = 81. From Fig. 7(a) to Fig. 7(c), the magnetic field in- 
creases gradually. We can see that in the weak magnetic field, 
the peaks of Sc are still regularly distributed for the zigzag 
ribbon and are irregular for the armchair ribbon due to the dif- 
ferent band structure for the zigzag edge and armchair edge 
ribbon. Moreover, for the insulating armchair edge ribbon, 
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the energy gap near Dirac point is diminished because of the 
magnetic field B, the very high and sharp S c at B = (see 
Fig. 5) is gradually dropped with the increasing of B. But at 
the weak magnetic field BSo/<f>o = 0.0005/7T, the Nc can still 
reach 3 (see Fig.7(a3)), which is much larger than all peaks 
of Sc in the high magnetic field case. Similar to Fig. 6, with 
the increasing of B further, the LLs is gradually formed from 
Dirac point to the high £>, the the properties of Sc for three 
systems gradually tend to the same. At the high magnetic field 
BSo/<f>o = 0.005/7T, LLs are completely formedfor \E F \ < 0.3, 
then Seebeck coefficient S c for three different systems are all 
the same to that in the Hall region. 

IV. CONCLUSION 

In summary, by using the Landauer-Biittiker formula com- 
bining with the non-equilibrium Green's function method, the 
Nernst effect in the crossed graphene ribbon and the Seebeck 
effect in the single graphene ribbon are investigated. It is 
found that due to the electron-hole symmetry, the Nernst coef- 
ficient Nc is an even function while the Seebeck coefficient S c 
is an odd function of the Fermi energy Ep. Nc and Sc show 
peaks when Ef crosses the Landau levels at high magnetic 
fields or crosses the transverse sub-bands at the zero magnetic 
field. In the strong magnetic field, due to the fact that high 
degenerated Landau levels dominate transport processes the 
Nernst and Seebeck coefficients are similar for different chi- 



rality ribbons. The peak height of Nc and Sc, respectively, 
are [ln2/|n|] and [ln2/«] with the peak number «, except for 
« = 0. For zeroth peak, it is abnormal. Its peak height is 
[2 In 2] for the Nernst effect and it disappears for the Seebeck 
effect. While in zero magnetic field, Nernst effect is absent 
and the Seebeck effect is strongly dependent on the chirality 
of the ribbon. For the zigzag edge ribbon, the peaks of S c are 
equidistance, but they are irregularly distributed for armchair 
edge ribbon. Surprisingly, for the insulating armchair edge 
ribbon, the Seebeck coefficient S c can be very large near the 
Dirac point due to the energy gap. When the magnetic field 
increases from zero to high values, the irregularly or regu- 
larly distributed peaks of S c in different chiral ribbons grad- 
ually tends to be the same. In addition, the nonzero values 
of the Nernst coefficient Nc appear first near the Dirac point 
and then gradually in the whole energy region. It is remark- 
able that for certain crossed ribbons, the Nernst coefficient Nc 
at weak magnetic fields can be much larger than that in the 
strong magnetic field due to small transmission coefficient in 
the transverse terminals. 



ACKNOWLEDGMENTS 

We gratefully acknowledge the financial support by a RGC 
grant (HKU 704308P) from the Government of HKSAR and 
NSF-China under Grants Nos. 10525418, 10734110, and 
10821403. 



1 K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, Y. Zhang, 
S.V. Dubonos, I.V. Grigorieva, and A.A. Firsov, Science 306, 666 

(2004) . 

2 K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, M.I. Kat- 
snelson, I.V. Grigorieva, S.V. Dubonos, and A.A. Firsov, Nature 
(London) 438, 197 (2005); Y. Zhang, Y.-W. Tan, H.L. Stormer, 
and P. Kim, Nature (London) 438, 201 (2005). 

3 C.W.J. Beenakker, Rev. Mod. Phys. 80, 1337 (2008); A.H. Castro 
Neto, F.Guinea ,N.M.R. Peres ,K.S. Novselov, and A.K. Geim, 
Rev. Mod. Phys. 81, 109 (2009). 

4 J. R. Williams, L. Dicarlo and C. M. Marcus, Science 317, 638 
(2008); A. f. Young and P. Kim, Nature Phys. 5, 222 (2009). 

5 K. Nakada, M. Fujita, G. Dresselhaus and M. S. Dresselhaus, 
Phys. Rev. B 54, 17954 (1996). 

6 K. Wakabayashi, M. Fujita, H. Ajiki and M. Sigrist, Phys. Rev. B 
59, 8271 (1999); 

7 V. P. Gusyninl, and S. G. Sharapov, Phys. Rev. B 73, 245411 
(2006). 

8 A. A. Abrikosov, Fundamentals of the Theory of Metals (North- 
Holland, Amsterdam, 1988); J. M. Iiman, Electrons and Phonons 
(Oxford University Press, Oxford, UK, 1960); C. W. J. Beenakker 
and A. A. M. Staring, Phys. Rev. B 46, 9667 (1992). 

9 M. Cutler and N. F. Mott, Phys. Rev. 181, 1336 (1969); D. K. C. 
Macdonald, Thermoelectricity (Dover, New York, 2006). 

A. S. Dzurak, C. G. Smith, C. H. W. Barnes, M. Pepper, L. Martin- 
Moreno, C. T. Liang, D. A. Ritchie, and G. A. C. Jones, Phys. Rev. 
B 55, R10197 (1997); R. Scheibner, H. Buhmann, D. Reuter, M. 
N. Kiselev, and L.W. Molenkamp, Phys. Rev. Lett. 95, 176602 

(2005) ; R. Scheibner, E. G. Novik, T. Borzenko, M. Konig, D. 



Reuter, A. D. Wieck, H. Buhmann, and L. W. Molenkamp, Phys. 
Rev. B 75, 041301 (2007). 

11 H.-C. Ri, R. Gross, F. Gollnik, A. Beck, R. P. Huebener, P. Wagner 
and H. Adrian, Phys. Rev. B 50, 3312 (1994); C. Hohn, M. Galffy 
and A. Freimuth, Phys. Rev. B 50, 15875 (1994); Z. A. Xu, J. Q. 
Shen, S. R. Zhao, Y. J. Zhang and C. K. Ong, Phys. Rev. B 72, 
144527 (2005); 

12 S. Onoda, N. Sugimoto, and N. Nagaosa, Phys. Rev. B 77, 165103 
(2008); Wei-Li Lee, S. Watauchi, V. L. Miller, R. J. Cava, and N. 
P. Ong, Phys. Rev. Lett. 93, 226601 (2004); T. Miyasato, N. Abe, 
T. Fujii, A. Asamitsu, S. Onoda, Y. Onose, N. Nagaosa, and Y. 
Tokura, Phys. Rev. Lett. 99, 086602 (2007). 

13 A. Banerjee, B. Fauque, K. Izawa, A. Miyake, I. Sheikin, J. 
Flouquet, B. Lenoir and K. Behnia, Phys. Rev. B 78, 161103(R) 

(2008) . 

14 N. Hanasaki, K. Sano, Y. Onose, T. Ohtsuka, S. Iguchi, I. 
Kezsmarki, S. Miyasaka, S. Onoda, N. Nagaosa and Y. Tokura, 
Phys. Rev. Lett. 100, 106601 (2008). 

15 K. Behnia, M. A. Measson, and Y. Kopelevich, Phys. Rev. Lett. 
98, 166602 (2007); K. Behnia, L. Balicas, Y. Kopelevich, Science 
317, 1729 (2007). 

16 Joshua P. Small, Kerstin M. Perez, and Philip Kim, Phys. Rev. 
Lett. 91, 256801 (2003). 

17 Y. M. Zuev, W. Chang, and P. Kim, Phys. Rev. Lett. 102, 096807 

(2009) . 

18 P. Wei, W. Bao, Y. Pu, C. N. Lau, and J. Shi, Phys. Rev. Lett. 102, 
166808 (2009). 

19 M. Cutler and N. F. Mott, Phys. Rev. 181, 1336 (1969); D. K. C. 
Macdonald, Thermoelectricity (Dover, New York, 2006). 



8 



N. M. R. Peres, J. M. B. Lopes dos Santos and T. Stauber, Phys. 
Rev. B 76, 073412 (2007). 

B. Dora and P. Thalmeier, Phys. Rev. B 76, 035402 (2007). 

H. Nakamura, N. Hatano, and R. Shirasaki, Solid State Commun. 

135, 510 (2005). 

R. Shirasaki, H. Nakamura, and N. Hatano, e-J. Surf. Sci. Nan- 
otechnol. 3, 518 (2005). 

M. Jonson and S. M. Girvin, Phys. Rev. B 29, 1939 (1984); H. 
Oji and P. Streda, Phys. Rev. B 31, 7291 (1985); J. Phys. C: Solid 
State Phys., 17, 3059 (1984)-3066. 

D. N. Sheng, L. Sheng, and Z. Y. Weng, Phys. Rev. B 73, 233406 
(2006); Z. Qiao and J. Wang, Nanotechnology 18, 435402 (2007). 
W. Long, Q.-F. Sun, and J. Wang, Phys. Rev. Lett. 101, 166806 
(2008); J. Li and S.-Q. Shen, Phys. Rev. B 78, 205308 (2008). 
In the references, [M. Jonson and S. M. Girvin, Phys. Rev. B 29, 
1939 (1984) and H. Oji and P. Streda, Phys. Rev. B 31, 7291 
(1985)], they pointed out that the standard Kubo formula can't 
apply to the Nernst coefficient in the presence of strong magnetic 
fields, because that a diamagnetic surface currents and a "diather- 
mal" surface currents appear under the strong mangetic fields. So 
a fundamental correction on the standard Kubo formula is neces- 
sary. However, the Landauer-Buttiker formula is to calculate the 
current flowing from a terminal to the center region, in which the 
both interior and surface currents have been included. So it is va- 
lidity regardless of the magnetic field. 

D. H. Lee and J. D. Joannopoulos, Phys. Rev. B 23, 4997 (1981); 
ibid, 23, 4988 (1981). 

S.-G Cheng, Y. Xing, Q.-F Sun and X. C. Xie, Phys. Rev. B 78, 
045302 (2008). 



P 

IV 



u 
o 
n 
a 
a 
to 



P 

c 



a: 
o 
o 

a 

0) 

« 



P 
o 



1.0 
0.5 
0.0 
-0.5 
-1.0 

0.5 

0.0 

-0.5 

0.5 
0.0 
-0.5 



(a1) 

I. 



(b1) 



0.002/jt 



; 0.001 1 

| o.005t 

V.'; — o.on 

w 



-+- 



-+- 



(d) 



0.005/Tt 



-o.ooit 

-•0.005t 
— 0.01t 
H 



V i.' 



001t 
005t 
01t 




-0.2 



0.0 



0.2 



Fermi Energy (t) 



-0.2 0.0 0.2 
Fermi Energy (t) 



-0.2 0. 
Fermi Er 



FIG. 7: (Color online) Panel (a)-(c) plot the Nemst coefficient Nc vs 
Fermi Energy E F at different temperatures k B T and different mag- 
netic fields BSa/<t>Q. The other parameters and the chirality of ribbon 
for the first, second, and third column panels are the same as Fig.4(a), 
(b), and (c), respectively. 



